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Context. A large fraction of stars host one or multiple close-in super-Earth planets. There is an active debate about whether these 
planets formed in situ or at greater distances from the central star and migrated to their current position. It has been shown that part of 
their observed properties (e.g., eccentricity distribution) can be reproduced by A-body simulations of in situ formation starting with a 
population of protoplanets of high masses and neglecting the effects of the disk gas. 

Aims. We plan to reassess the in situ formation of close-in super-Earths through more complete simulations. 

Methods. We performed A-body simulations of a population of small planetary embryos and planetesimals that include the effects 
of disk-planet interactions (e.g., eccentricity damping, type 1 migration). In addition, we also consider the accretion of a primitive 
atmosphere from a protoplanetary disk. 

Results. We find that planetary embryos grow very quickly well before the gas dispersal, and thus undergo rapid inward migra¬ 
tion, which means that one cannot neglect the effects of a gas disk when considering the in-situ formation of close-in super-Earths. 
Owing to their rapid inward migration, super-Earths reach a compact configuration near the disk’s inner edge whose distribution of 
orbital parameters matches the observed close-in super-Earths population poorly. On the other hand, simulations including eccentric¬ 
ity damping, but no type I migration, reproduce the observed distributions better. Including the accretion of an atmosphere does not 
help reproduce the bulk architecture of observations. Interestingly, we find that the massive embryos can migrate inside the disk edge 
while capturing only a moderately massive hydrogen/helium atmosphere. By this process they avoid becoming giant planets. 
Conclusions. The bulk of close-in super-Earths cannot form in situ, unless type I migration is suppressed in the entire disk inside 1 
AU. 

Key words. Planets and satellites: formation - Planets and satellites: atmospheres - Planet-disk interactions - Methods: numerical 


1. Introduction 

Recent observations of exoplanets have revealed a large num- 
O ber of close-in lo w-mass planets (e.g., ISchneider et al.l |201 lb 
IWright et al.l201 ih . As of January 2015, 337 systems harbor 839 
OQ planets with masses M < lOOMe (or with radii R < lOR©) and 
semimajor axes a < lAU (or with orbital periods P < 200day). 
We define these as “close-in super-Earths.” 

^ They have a semimajor-axis distribution centered on 0.1 AU. 
Super-Earths in each system are generally confined within a 
^ few tenths of an AU. Although some of these close-in super- 
. . Earths are in or near first-order mean motion resonances (es- 
pecially in 3:2 resonances), a large number of cl ose-in super- 
Earths are not in rnean m otion resonances (e.g., iMavor et al.l 
l2009tlLrssauer et al.l201 lb|) . Typical orbital separations between 
a these planets are 10 - 30rH, where rn is the mutual Hill radius, 
which is similar to orbital separations between the solar system 
terrestrial planets. Period ratios of adjacent pairs lie between the 
4:3 and 3:1 resonances, respectively. It is estimated that a large 
number of super-Earths have small eccentricities e ~ 0.01 - 0.1, 
while some of them can hav e high eccentricities up to e ~ 0.5 
(e.g.. iMoorhead et al.l 1201 ll) . The mutual inclinations between 
planetary orbits could be estimated in a fraction of transiting 
systems from the K epler catalog and app ear to be low with an 
average of i ~ 0.03 dFabrvckv et al.ll2()3 . 

iHansen & Murravl (l2012Ll2013l) present A-body simulations 
of the in situ formation of close-in super-Earths from embryos 
that are placed between 0.05 AU and 1 AU. To account for 
the masses of the known planetary systems, they had to assume 


that up to 100 Earth masses of solids existed in the disk within 
1 AU, implying a surface density of solids much higher than 
for the minimum mas s solar nebula (MMSN) (IWeidenschillingl 
I 1977 I : iHavashil Il98lh . Assuming a nominal gas/solid density 
ratio of 100, this would imply a very massive protoplanetary 
disk that would probably be g ravitationally unstable. However, 
IChatteriee & TanI (l20l4U()T^ propose that the disk may be en¬ 
riched in solids relative to the gas thanks to the inward migration 
of dust grains, pebbles, and planetesimals. Thus, it may be legit¬ 
imate to assume a very high surface density of solids embedded 
in the inner protoplanetary disk with a mass of gas comparable 
to gas in the MMSN model. 

The simulations of IHansen & Murravl (1201 2i 12013h . never¬ 
theless, are quite simplistic. They start from a population of 
protoplanets of high masses (sometimes multiple Earth masses), 
which are supposed to have alr eady achieved the com pletion of 
the oligarchic growth process (iKokubo & Idalll998l) . No plan¬ 
etesimals are considered, and the effects of the gas, in relation to 
the migration and eccentricity or to the inclination damping of 
the protoplanets’ orbits, are not taken into account. With these 
assumptions, their simulations suggest that an in situ accretion 
without gas can explain the distributions of orbital periods and 
eccentricities of observed super-Earths. However, their simula¬ 
tions did not reproduce the distribution of period ratios between 
adjacent super-Earths, because t hey had a deficit of close pairs 
and resonant pairs (see Eig. 15 in IHansen & Murr^ (120131) ). 

Int erestingly, the same deficit was found by ICossou et al.l 
(I 2 OI 4 I) in a different model accounting for the migration of plan¬ 
etary embryos from several AUs away. On the other hand, a 
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model by lOgihara & Idal (l2009h with gas drag and type I mi¬ 
gration led to more separated non-resonant pairs, when type I 
migration was reduced by about a factor of 100. 

In this paper, we revisit t he process of in situ formation 
(iHansen & Murravl2012ll^l3h with more complete simulations 
that start from a population of small planetary embryos (i.e., 
Mars-mass) and planetesimals, both carrying cumulatively 50% 
of the solid mass in the disk. These initi al conditions are typ - 
ical of terrestrial-planet simulations (e.g., lO’Brien et al]l200^ . 
Thus, we do not assume that the protoplanets have reached 
the completion of the oligarchic growth process, but we sim¬ 
ulate that process from a much more primordial state. More¬ 
over, we consider the action of the gas that forces a migration 
and eccentricity/inclination damping of the embryos and plan- 
et esimals. Our aim is to c l arify differences from the studies 
of iHansen & Murr^ (1201 2i 1201 3h and to understand, with im¬ 
proved and more realistic simulations, whether in situ formation 
can explain the observed properties of the close-in super-Earth 
systems. We performed simulations with different disk models 
(i.e., the amount of solids available), but we also present a model 
without migration and one without gas for comparison. 

In addition, we extend our model by including the accre¬ 
tion of primitive atmospheres onto super-Earths. The analyses of 
transit observations with radial velocity measurements or analy¬ 
ses of transit-time variations show that most of the planets larger 
than 2.5 Earth radii have v ery low densities (iMarcv et al.ll20l4 
iHadden & Lithwi^l2014ll . so are thought to have thick H/He 
atmospheres (up to 10-20% by mass). We therefore consider 
the acquisition of H/He envelopes and discuss whether observed 
low-density super-Earths can be explained by our in situ accre¬ 
tion model. 

The rest of the paper is organized as follows. In Sect. 
we describe our model and methods of A^-body simulations; in 
Sect. [3 we give a series of the results of our simulations; in 
Sect. ID we present results of simulations that include accretion 
of H/He atmosphere; in Sect.|5]the conclusions are provided. 

2. Model and methods 

2.1. Disk model 

We assume that the gas distribution is similar to the one of the 
classical MMSN model. Thus, for the gas surface density, we 
assume 

where 2g and r are the gas surface density and the radial distance 
from the central star, respectively. The value of Eg is 1.4 times 
the value in the MMSN. The gas dissipation is modeled as an 
exponential decay with the depletion timescale of 1 Myr. We set 
the disk inner edge at r = 0.1 AU. The inner edge of the disk is 
expected to be at the radius where the orbital period equals the 
stellar rotation period. This could be a factor of two-three smaller 
than we assume. Our assumption of an inner edge at 0.1 AU is 
dictated by constraints due to computational time. Our results 
concerning the position of the final planets relative to the disk’s 
inner edge can scale with the assumed edge’s position. The tem¬ 
peratu re distribution is that for an optically-thin disk (iHavashil 
EMI, so that the disk scale height is 

/ r \ 5/4 


Table 1. List of models. In Model 2, the initial solid amount is decreased 
by a factor of two from Model 1. Model 4 corresponds to the model of 
IHansen & Mu^ 


Model 

Type I migration 

e, /-damping 

1 

yes 

yes 

2 

yes 

yes 

3 

no 

yes 

4 

no 

no 


2.2. Initial conditions and numerical method 

We chose an in i tial s olid distribution similar to that of 
IHansen & Murravl (l2012l) : that is, 50M^ in total are placed be¬ 
tween 0.1 and lAU. In our standard model, we set 250 embryos 
with a mass of M - O.IM 0 and 1250 planetesimals with a mass 
of M - O.O 2 M 0 in such a way as to keep the radial distribution 
of the solid surface density proportional to Planetesimals 
gravitationally interact with the embryos but not with each other. 
This set-up is typical of successful simulations for the growth of 
the terrestrial planets in our solar system. By adopting the same 
set-up, we explore the effects of the enhanced solid distribution 
and much shorter orbital distances. 

There may be some caveats associated with the initial condi¬ 
tion. As shown in Sect. [2 once a disk of planetesimals and plane¬ 
tary embryos is set, the growth of planets in the close-in region is 
quite rapid. If the planetesimal formation process takes a longer 
time than the planet-growth timescale, it would be important to 
take planetesimal formation into account in simulations. How¬ 
ever, the formation of planetesimals is not yet fully understood. 
Therefore, in this study, we do not attempt to include the plan¬ 
etesimal formation; instead, we start simulations with already 
formed planetesimals. We note, however, that it is unlikely that 
planetesimal formation takes a timescale comparable to the disk 
lifetime (~ Myr) for two reasons. Eirst, this timescale would cor¬ 
respond to 30 million orbital periods, and it is difficult to under¬ 
stand why it should take so long. Eor instance, in the solar sys¬ 
tem chondritic planetesimals formed in about one million orbits. 
Second, the planetesimal formation process is mo st likely related 
to the drift of small particles through the disk (iJohansen et al.l 
l2014l) . which is due to gas drag, so very likely most of the mass 
was fed to the inner disk at an early time, when the gas density 
was higher. A late formation of the planetesimals would require 
that the small particles drift into the inner part of the disk only 
when the disk is disappearing. We think that this possibility is 
difficult to envision. 

Table [1] shows the list of simulations for each model. In 
Model 2, the total mass inside lAU is decreased by a factor of 
two (the total mass is 25M®). In Model 3, type I migration is ne¬ 
glected, but eccentricity damping is still considered. In Model 4, 
the effect of gas is ignored th e same as in the simulation by 
IHansen & Mum\^(l2012ll2013h . 

Our A-body code is based on SyMBA (Duncan et al. 1998), 
modified so that the effects of the gas disk are included accord¬ 
ing to the formulae reported in Sect. 12.31 When bodies collide 
with each other, the bodies are merged by assuming perfect ac¬ 
cretion. The physical radius of a body is determined by its mass, 
assuming an internal density of p = 3 g cm“^. The inner bound¬ 
ary of the simulation is set to r = 0.05AU. We use a 0.0004-year 
timestep for integrations. 
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2.3. Effects of gas 

Eccentricities, inclinations, and semimajor axes are damped by 
disk interaction. Planets with more than roughly 0.1 Me suffer the 
tidal damping by the density wave, while planetesimals undergo 
aerodynamical gas drag. 


2.3.1. Damping for embryos 

The eccentricity damp ing timescale for embryos, 4 , is given by 
dTanaka & Wardl2004 


4 = 
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where Cd, fp,Pg, and Ah are the gas drag coefficient, the physical 
radius, the density of gas disk, and the velocity of the body rela¬ 
tive to the disk of gas, respectively. Although our planetesimals 
have a mass of 0.02 Earth masses, we consider them as super- 
planetesimals, representing a swarm of much smaller objects cu¬ 
mulatively carrying the same total mass. Erom the size distribu¬ 
tion of the asteroid belt, we thin k that most planetesim als had a 
physical size of 50km in radius (iMorbidelli et al.ll2009l) . Eor the 
calculation of gas drag, this is the size we assume. Eor the value 
of Cn, we use the same definition as in previous studies (e.g., 
lAdachi et al.lll97^ Brasser et al. 2007), which depends on the 
Mach number, the Knudsen number, and the Reynolds number. 


3. Results 

3.1. Outcomes of Models 1 and 2 


where M, is the stellar mass, L* the stellar luminosity, Cs the 
sound speed, uk the Keplerian velocity, and Q. the orbital fre¬ 
quency, respectively. Here the relative motion between gas and 
planets is assumed to be subsonic (evj^ < Cs). Eor planets with 
high eccentricities and inclinations, we include a correction fac¬ 
tor according to Eqs. (11) and (12) of Cresw ell & Nelson (2008) . 

The migration timesc ale 4 is given by (iTanaka et al.ll200a 
iPaardekooner et aTll201 ih 



where jS is a coefficient that determines the direction and speed 
of type I migration. The type I migration torque depends on the 
Lindblad torque, the barotropic part of the horseshoe drag, the 
entropy-related part of the horseshoe drag, the barotropic part 
of the linear corotation torq ue, and the entropy - related part of 
the linear corotation torque. IPaardekooner et aP (1201 ll) derived 
the total type I migration torque, including both saturation and 
the cutoff at high viscosity. We write the migration coefficient fi 
entering in Eq. (HI in the form 

(3 —^c.baro (5) 

where /^l, jSc.baio, and/3c,ent are related to the Lindblad torque, the 
barotropic part of the corotation torque, and the entropy-related 
part of the corotatio n torque, respectively . Each formula is given 
by Eqs. (11)-(13) in lQgitiara et al.l (l2015h . In addition, the coro- 
tation torque decreas es as the planet eccentricity increases (e.g., 
iBitsch & Klevll2010ll. We also consider th is effect using the fol¬ 
lowing formul ^Eendvke & Nelsonll2014l) : 

y6c,baioexp|-^j, (6) 

ySc,entexp|-^j, (7) 

where Cf = 0.5/r/r + 0.01. 


^C,baro(^) “ 
^C,ent(^) ~ 


2.3.2. Damping for planetesimals 

The a erodynamical gas drag force per unit mass is (lAdachi et al.l 
119761) 


F 


aero 


-^CoTTr^PgAMAH 


( 8 ) 


The formation of close-in super-Earth systems can be schemati¬ 
cally divided into three phases: (i) growth of embryos from plan- 
etesmals, (ii) migration of embryos, and (iii) gas depletion and 
long-term orbital evolution. Through a series of simulations, we 
find that the first accretion stage is extremely rapid (~0.1 Myr, 
much faster than the typical timescale of several 10 Myr ob¬ 
served in terrestrial planet simulations). Thus, the effect of the 
enhanced amount of solid is not just that of producing bigger 
planets: it also reduces the formation timescale. The closer prox¬ 
imity to the central star (i.e., sh orter orbital peri ods) also favors a 
much faster accretion (see also lLee et al.ll2014l) . Eigure[T] shows 
snapshots of the evolution of one simulation for Model land 
indicates the gas surface density at each time (right axis). Eig- 
ure| 2 |a) shows the time evolution of the semimajor axis for this 
run. The color of lines indicates the eccentricity of the planets 
(see color bar). At f = lO^yr, almost all planetesimals initially 
placed inside r ^ 0.2AU have been accreted by embryos. The 
first accretion phase ends before t - O.OlMyr and O.lMyr inside 
r ^ 0.4AU and ^ 0.7AU, respectively. 

Because of the short accretion timescale, it is thus not 
correct to neglect the gas effects as in Hansen and Murray’s 
works. In fact, the protoplanets become massive well before 
the gas disk is substantially depleted. The gas forces the plan¬ 
ets to migrate inward. All planets would be lost into the 
star if there were no inner edge of the disk. With a sharp 
disk inner edge, the inn ermost planet is tra pped at the edge 
by the planet-trap effect dMasset et al.l 1200^ . The other plan¬ 
ets pile up in mutual mean motion resonances with the for¬ 
mer. The typical commensurabilities are between 5:4 and 6:5 
with orbital separations of ^ 5 - lOre, thus the final orbits 
are packed near the disk’s edge. The resonant configurations 
of two bodies that undergo convergent migration are deter¬ 


mined by the mass of the bodies and the relative miratio n 
speed (e.g..lMustill & Wvatl 1201 lUOgihara & Kobavashil2013l) . 


lOgihara & Kobavashil ( 2013 ) have derived the critical migration 
timescale for capture into first-order mean motion resonances. 
They found that if the relative migration timescale is shorter than 
fa.crit - 1 X 1 O^(Mi/M 0 )“"^^^ 7 ’k, where M\ is the mass of larger 
body and Tk the Keplerian period, bodies can only be captured 
in reson ances closer than the 4:3 reso nance (see Eq. ( 6 ) and Ta¬ 
ble 2 in lOgihara & Kobavashill2()T^ . In the result of Model 1, 
the mass of migrating embryos is ^ O. 5 M 0 and the mass of plan¬ 
ets near the edge is ~ IM 0 , so the relative migration timescale 
is - IO^Tk (Eq. (H). Migrating embryos therefore settle into 
closely packed configurations. Live planets form at the end of the 
simulation presented in Eigs. [T]and|2|a). In this run, the planets 
do not exhibit orbital instability after gas depletion because the 
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Fig. 1. Snapshots of a system for Model 1. Filled circles represent bod¬ 
ies. The size of the circles is proportional to the radius of the body. The 
smallest circle represents a 0.02 Earth-mass body, while the largest one 
represents a 33 Earth-mass body. The solid line indicates the gas surface 
density (right axis). 


number of planets in the system is small and the planets are in 
me an motion resonances, leading to a long-lasting orbital stabil¬ 
ity ([Chambers et al.lflO^lMatsumoto et alj|2012l) . 

We performed ten simulation runs for each model with dif¬ 
ferent initial positions of embryos and planetesimals. The results 
are qualitatively the same: final orbits are compact near the disk’s 
inner edge. In some runs, the system undergoes orbital instability 
during the third phase after 1 Myr, resulting in non-resonant and 
relatively separated orbits with a smaller number of planets (see 
Fig. |2jb) for example). Figure |3l a) shows the final orbital con¬ 
figurations of all ten runs, where the planets that formed through 
the same run are connected with a line. We observe steep mass 
gradients in this figure; the largest bodies are located near the 
edge, and the planetary mass monotonically decreases when the 
semimajor axis increases. This is because, in the presence of a 
strong migration torque, the resonant system is stable only if 
the innermost planet is the most massive, as already found by 
iMorbidelli et al.l (l2008h . If originally the innermost planet is not 
the most massive, an instability typically occurs. The first and 
the second planets have encounters with each other, and the sys¬ 
tem stabilizes in resonance only after that the two planets have 
exchanged their relative positions. The same is true for the sec¬ 
ond, relative to the third planet and so forth. 

Figure [3b) shows orbital configurations of observed close- 
in super-Earth systems, in which the mass and semimajor axis 
of all planets are known. It is clear that the bulk architecture 
of the observed systems is inconsistent with the simulated plan¬ 



Time (yr) 


Fig. 2. Time evolution of planets for Model 1. The filled circles con¬ 
nected with solid lines represent the sizes of planets. The smallest circle 
represents a 0.2 Earth-mass embryo, while the largest ones represent a 
33 Earth-mass planet in panel (a) and 35 Earth-mass planet in panel (b). 
The color of line indicates the eccentricity (color bar). 


etary systems. No steep mass gradients are observecQ, and the 
bulk architecture of the observed systems cannot be reproduced 
through the simulations of model 1. If embryos migrate from 
beyond I AU, the mass gradient would be shallower and/or the 
outer boundary of the final planet distribution at around 0.2 AU 
would be removed. This would suggest that a “migration model” 
yields better results in reproducing the observed close-in super- 
Earth systems, which should be investigated by future A-body 
simulations. 

Eigure|4|a) shows period ratios of adjacent pairs of observed 
close-in super-Earths indicating period ratios of first-order mean 
motion resonances (e.g., 2:1 and 3:2) and the cumulative dis¬ 
tribution of the period ratios (right axis). Most pairs have period 
ratios between 4:3 and 3:1, and there are a few pairs that have pe¬ 
riod ratios lower than 4:3. Eigurellfb) shows the period ratio dis¬ 
tribution for the ten runs for Model 1 and a copy of the observed 

' The mass gradient is steep in the Kepler-101 system, where the in¬ 
nermost planet has 51 Earth mass and the second planet is about 3 Earth 
mass. 
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Period Ratio 


Fig. 4. Comparison of the distributions of period ratios of adjacent pairs of planets for (a) observation and (b)-(d) simulations. The distribution 
of period ratios is presented as a histogram (see the left y-axis) and as a cumulative distribution (see the right y-axis). Panels (b), (c), (d), and (e) 
show the results of Models 1, 2, 3, and 4, respectively. The dashed lines in panels (b)-(d) represent the cumulative distribution of observed planets 
shown in red in panel (a). The vertical lines indicate locations of mean motion resonances. 


distribution. Some pairs are in closely-spaced resonances (e.g., 
5:4), and others have been knocked out of resonant orbits during 
the late instability. Although very closely spaced pairs (< 4:3) 
can form in Model 1, the cumulative distribution of observed 
close-in super-Earths is not matched by the results of Model 1. 
In fact, a Kolmogorov-Smirnov (K-S) test indicates that the cu¬ 
mulative distributions are statistically different (2ks ^ 0.01). 

Figure |5] shows cumulative distributions of the eccentricities 
of the observed close-in super-Earths and of all the results pro¬ 
duced in the simulations corresponding to the same model. The 
general trend is that the eccentricity of the results of Model 1 is 
smaller than for close-in super-Earths. In Model 1, planets with 
e < 0.03 account for 66 percent of all bodies, while exoplanets 
with e < 0.03 make up only 26 percent of all planets. One reason 
for the small eccentricity in Model 1 is that the number of plan¬ 
ets in a system is small (A^ - 3 - 6), and the orbital stability time 
is long (IChambers et al.l[T996h . which inhibits giant impacts dur¬ 
ing gas dispersal (see Fig.|2la) for example). In addition, even if 
planets undergo orbital instability as in Fig.|2b), the eccentricity 
is not highly excited. This is because the effect of mutual scat¬ 
tering between planets is limited by the small number of planets. 
Eccentricity damping also operates during the gas dissipation 


phase due to remnant gas (see Fig. |2|b), for example). A K-S 
test indicates that the observed and simulated eccentricity distri¬ 
butions are statistically different (Qks ^ 0.01). In summary, the 
results of Model 1 cannot reproduce the bulk properties (period 
ratio, eccentricity) of observed close-in super-Earth systems. Ec¬ 
centricities of exoplanets can be overestimated (see Sect. l3.2l for 
discussion). Thus, our results disagree with those of Hansen and 
Murray, which is not surprising given that the latter neglected 
the effect exerted by the disk of gas (particularly the inward mi¬ 
gration of proroplanets), not realizing that the super-Earths must 
form well within the gas-disk lifetime. 

We also performed ten runs for Model 2, where the initial 
solid amount is reduced by a factor of two from Model 1, in 
the hope of observing a slower accretion rate and consequently 
weaker migration effects. However, the results are qualitatively 
the same with those of Model 1. In Fig.|4jc), planetary pairs with 
their period ratio of < 4:3 account for 78 percent of all pairs, 
and the cumulative distribution is different from that of observed 
close-in super-Earths. The results are even more closely packed 
than those of Model 1 because planets are less vulnerable to or¬ 
bital instability during the third phase. The eccentricity distribu¬ 
tion in Fig.|5]also differs from that of close-in super-Earths. 
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Semimajor Axis (AU) 



Semimajor Axis (AU) 


Fig. 3. Results of 10 simulations of Model 1 (panel (a)). Observed close- 
in super-Earths systems (panel (b)). 



Time (yr) 

Fig. 6. Same as Fig. |2] but for a representative simulation of Model 3 
(panel (a)) and Model (4) (panel (b)). 


3.2. Outcomes of Models 3 and 4 



Fig. 5. Comparison of cumulative eccentricity distributions between ob¬ 
served close-in super-Earths (thin solid line) and the planets formed 
through simulations (thick lines, see legend). The thin dotted line indi¬ 
cates the cumulative eccentricity distribution of observed super-Earths, 
in which each eccentricity is assumed to be e - cr. 


We then present results of Models 3 and 4, in which we sup¬ 
press the migration torques exerted by the gas-disk interaction 
onto the planets (Model 3) or we neglect the presence of gas 
altogether (Model 4). Clearly, both models are academic. In 
fact, a general mechanism that suppresses type I migration has 
never b een found. Locally, type I migra t ion can be halted o r re¬ 
versed (iPaardekooner & MellemalbOOhl : iBitsch et'^l2014l) . but 
no global weakening of type I migration has even been demon¬ 
strated. As for the absence of gas, this seems inconsistent with 
the fast growth timescale for the super-Earths. It is difficult to 
imagine that the gas disappears signihcantly faster than what 
we assume above (1 Myr). The reason we present these models 
is to highlight the role of migration and eccentricity/inclination 
damping in the results we presented before. 

Figures |6ja) and (b) show the typical orbital evolution for 
Models 3 and 4, respectively. In both cases, planetary systems 
are not as compact when compared with the results of Models 1 
and 2. In the simulations for Model 3, in which type I migra¬ 
tion is neglected, planets undergo slow in ward migration due t o 
eccentricity damping (see also Sect. 5.2 in lOgihara et al.ll2014l) . 
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and planets are temporarily captured in mutual mean motion res¬ 
onances before a few Myr. Then, they undergo orbital instabil¬ 
ity and collide with each other, resulting in non-resonant orbital 
configurations, which are qualitati vely similar to those o btained 
in the slow-migration simulation of lQgihara & Idal(l2009l) . In the 
results for Model 4, planets are never in resonances, which is 
alm ost the same as in the simulati ons of ICossou et al.l (l2014h 
and iHansen & Muira^ (l2012i l2013h . The averaged mass of the 
largest bodies is ^ 13 Me (Model 3) and ^ 14 Me (Model 4), 
which are lower than for Model 1. 

Figures |4jd) and (e) show the period ratio distribution of 
Models 3 and 4, respectively. Interestingly, the results match the 
observations much better than those of Models 1 and 2. In the 
results of Model 3, most pairs have period ratios between 4:3 
and 3:1, while some pairs are relatively separated (> 3:1). In the 
results of Model 4, 85 percent of pairs lies between 2:1 and 3:1. 
In comparison with the observed distribution. Model 3 is a good 
match to the distribution of exoplanets. A K-S test shows that the 
distributions are similar with a significance level of Qks = 0.24 
for Model 3, while the distribution of Model 4 is not similar to 
the observed distribution (2 ks ^ 0 .01). 

In Fig.|5] the eccentricity of the planets produced in Model 3 
is generally smaller than in Model 4 because of the eccentric¬ 
ity damping. Compared with the distribution of exoplanets, the 
K-S test shows that the distributions are different for Model 3 
(2 ks ^ 0.01), while the distributions for Model 4 are closer, but 
are still not satisfactory (2 ks = 0.026). Thus, neither Model 3 
nor Model 4 explains the observation, because the first has prob¬ 
lems with the eccentricity distribution, the second with the or¬ 
bital period distribution. 

However, it is fair to say that the measurement of the eccen¬ 
tricity of exoplanets is still difficult, and the uncertainty on the 
results is quite large. In particular, it has been shown that ec¬ 
centricities, which are d erived from radial velocity surveys, can 
be ov erestimated (e.g., IShen & Turned l2008t IZakamska et al.l 
1201 ih . Therefore, the observed eccentricity distribution in Fig.|5] 
may shift to lower values. In this case the results of Model 3 
might match the eccentricity distribution as well. 

As an example, we recalculated the eccentricity distribution 
in a simple way. The eccentricity of each exoplanet is set to e-cr, 
where cr is the estimated error. The new distribution is indicated 
in Fig. |5] We find that the new distribution gives a better match 
to Model 3 (Qks = 0.30) rather than Model 4 (2 ks = 0.010). 


4. Accretion of primitive atmospheres 


accretion to compute the mass of atmosphere expected for our 
super-Earths. 


4.1. Model 


Once planetary embryos embedded in a gas disk are sufficiently 
massive (typically the mass of Mars or more), they capture 
part of the disk gas to have an atmosphere of their own (e.g., 
IWuchterl et alj[2000l) . This atmosphere or envelope grows with 
the mass of the embryo itself until it becomes greater than the 
mass of the solid core. At this so-called crossover mass (Menv ~ 
Mcore), accretion enters a runaway phase with the envelope mass 
increasing exponent ially and the planet becoming a giant planet 
(iPollack et al.lll99^ . In this process, the accretion of solids has 
two effects: it increases the mass of the planet and heats the enve¬ 
lope. The first favors the growth of the envelope, but the second 
leads to an increase in the crossover mass because of a more ten¬ 
uous envelope. Most works on the growth of giant planet cores 
have thus focused on obtaining exp ressions that depen d on the 
accretion rate of planetesimals (e.g.. llkoma et al.lbOOlh . 

However, since in our simulations the accretion of solids 
stalls while the gas disk is still present, at this later stage, the rate 
of cooling of the envel ope becomes crucia l for controlling the 
growth of the envelope (iPollack et al.lll99^ . We choose to em- 
pi rically model the grow t h of the envelope by fi tting t he results 
of llkotna & Horil (1201 2l) . iPiso & Youdin (l2014h . and iLee et all 
(l2014h . which are works that account for the planetary cooling to 
calculate the resulting envelope growth. Specifically, we model 
the envelope mass, Menv, as 


Me, 


Me, 


ki 

l+k3 


1 +k3 



^^l(f/frun 1) 


(9) 


where Mcore is the core mass and the time to get to 
the crossover mass. As for the coefficients, we use ki = 
(Mcore/15M®)^', k 2 = I + (Mcore/40Me), and ki - 9. Most of 
the uncertainty is linked to the value of fnm, which depends crit¬ 
ically on the opacities chosen, the cooling rate of the core, and 
orbital distance. For the present simulations we choose to use as 
a fiducial value 


f„„ = 10"(Meore/5M®)-3 yr, (10) 

which lies in between the resu lts o btained by Ikoma & Horil 
(1201 2h and iPiso & Youdinl(l2014h and lLee et al.f(l2014h . 

Then the accretion rate onto the planet is expressed by 


An objection often expressed against the in situ accretion model 
for super-Earth s is that observations in dicate in many cases (e.g., 
Kepler-11, see iLissauer et al.ll201 lah that these planets have a 
low bulk density. One would expect planets grown in the inner 
part of the disk to be rocky, given that the high local disk tem¬ 
perature should not have allowed ice condensation. 

A possibility, however, is that super-Earths accreted prim¬ 
itive H/He atmospheres, leading to low bulk densities. Recent 
simulations of the structure and evolution of planetary atmo¬ 
spheres have demonstrated that super-Earths can indeed accrete 
primordial atmospheres from the protoplanet ary disk provide d 
that they do not accrete solids at a high rate (iLee et al.ll2014ll . 
As we have seen above, we expect that in situ formation is ex¬ 
tremely rapid and reaches completion well within the lifetime 
of the gas disk. Thus, we expect our planets to be in the condi¬ 
tion of a very low accretion rate of solids when there is still gas 
in the disk, enabling the acquisition of significant atmospheres. 
In what follows we implement the most recent recipes for gas 


Menv 


— Menv 


ki 1 

'^&+kit 



( 11 ) 


Equation (0 implicitly assumes that the disk can supply all the 
gas that the planet is able to accrete. In reality, this amount is 
limited by t he viscous inflow in the disk at the location of the 
planet (e.g.. lTanigawa & Ikomal2d(y^ . The accretion rate due to 
viscous diffusion is 


Mvis - 37rvZg, (12) 

where an “alpha model” for the disk viscosity is used with v = 
acji (Cs is the isothermal sound speed), and we adopt a = 10“^. 
The actual accretion rate is limited by the minimum of Eqs. (fTTT i 
and (fT2Ti . 

Based on these prescriptions, we recalculate our A-body sim¬ 
ulations for Model 1, this time accounting for the accretion of an 
atmosphere after the end of core accretion phase (t < lO^yr). 
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Fig. 7. Evolution of the semimajor axis and the envelope mass (color 
bar). Panel (a) shows the result of = lO’CMcore/SMe)^^ yr, while 
panel (b) indicates that of = lO'’(Mcore/5M0)“^ yr. The filled circles 
connected with solid lines represent the sizes of planets. The largest 
circles represent a 33 Earth-mass planet in panel (a) and 36 Earth-mass 
planet in panel (b). 


4.2. Results 

We now present the results of these simulations with the set-up 
of Model 1 but accounting for atmosphere accretion. There are 
two crucial questions that we wish to address; 1) Can atmosphere 
accretion be substantial enough to significantly reduce the appar¬ 
ent bulk density of the planet and mimic low-density objects as 
in the Kepler-11 system? 2) Can atmosphere accretion change 
the masses of the planets enough to induce a late dynamical in¬ 
stability that can result in less compact systems. 

Figure|2]shows the orbital evolution of each system in which 
the color of each planet’s lines indicates the ratio of enve¬ 
lope mass to total mass, Menv/ATtot (see color bar). We adopt 
bun = lO^jMcore/SMe)^^ yr as a fiducial value for trun in the 
result of panel (a). In addition, we also perform a simulation 
under more efficient conditions for envelope accretion fmn = 
lO®(Mcore/5M0)“^ yr, the results of which are shown in panel 
(b). 


In Fig.|2ta), the planet at a = 0.1 AU with Mcore - 9.2 Mq 
accretes gas from the disk and ends up retaining a thick atmo¬ 
sphere of Menv = 6.2 Mq and Menv/ATtot = 0.30 at f = lOMyr. 
This planet migrates inside of the disk inner edge at f ^ 2.6Myr 
by the interaction with four outer bodies. This migration prevents 
the former planet from accreting more gas. The innermost planet 
is also moved inside the disk inner edge before t - O.lMyr and 
stops accreting gas at that point. 

In the simulation with a shorter value of shown in 
Fig- mb), a planet with Menv = 14M® and Menv/Mtot = 0.56 
eventually forms and also moves inside of the disk edge at 
t ^ 0.3Myr. It is interesting to notice that as the envelope mass 
increases, systems become vulnerable to orbital instability. In 
panel (b), in fact, the system undergoes orbital instability earlier 
than in panel (a), and the number of final planets is also smaller. 

As for the first question posed in the beginning of this sec¬ 
tion, we find that in both simulations, one planet can acquire a 
thick H/He atmosphere from the disk, which may explain the ori¬ 
gin of the observed low density super-Earths. Regarding the sec¬ 
ond question, we also observe that the acquisition of a massive 
atmosphere by one planet destabilizes closely spaced systems, 
leading to relatively well separated systems with fewer planets 
or even single-planet systems. 

The properties (e.g., orbital separation and bulk density) of 
the final systems shown in Fig. |7Ja) are reminiscent of some 
known planetary system. For example, two super-Earths were 
discovered in the Kepler-36 system, where the inner planet 
would be a rocky planet without a thick atmosphere, and the 
outer one would possess a thick atmosphere. However, the typi¬ 
cal properties of the observed close-in super-Earths are unlikely 
to be reproduced. This is because the systems are tightly packed 
near the edge before the acquisition of atmospheres (see f = 10^ 
yr in Eig.|7]l, so that systems separated more than 2:1 resonances, 
which have been observed in super-Earth systems (see Eig.|4ja)), 
can hardly be reprod uced by simulations for Model 1. 

iLee et al.1 (1201 4l) point out that super-Earths with a mass of 
IOM 0 tend to undergo runaway gas accretion, thus becoming gas 
giant planets. We do not observe this phenomenon in the results 
of Eigs. |7ja) and (b), even though the core mass is high. This 
is because the massive planets move inside the disk inner edge, 
where there is no gas, and cease envelope accretion. 


5. Discussion and conciusions 

We have re-examined in situ formation of close-in super-Earths 
with improved simulations, in which the effects of the disk of 
gas are considered. The simulations are started with small em¬ 
bryos and planetesimals. We find that the accretion of planets 
is extremely rapid owing to the large amount of solid material 
and short orbital periods. Thus, it is not correct to ignore the ef¬ 
fects of the gas disk for the investigation of in situ formation of 
close-in super-Earths. 

We performed ten runs of simulation for our fiducial model 
and found the following. 1) The orbital architecture of resultant 
systems is very compact near the disk inner edge. 2) The eccen¬ 
tricity of super-Earths is small because planets can be stable after 
gas depletion. If they undergo orbital instability, the eccentricity 
is not highly excited. 3) The masses of planets monotonically 
decrease when increasing the semimajor axis. These character¬ 
istics are not consistent with observed close-in super-Earths. In 
fact, the cumulative observed distributions of period ratios of ad¬ 
jacent pairs and of eccentricities are statistically different from 
those we produce. 
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We have also investigated orbital evolution including the ac¬ 
cretion process of primitive atmospheres onto the super-Earths. 
The results show that close-in super-Earths that formed in situ 
can acquire a thick H/He atmosphere in which the planets stop 
envelope accretion when they migrate inside the disk’s inner 
edge. 

Interestingly, if type I migration is neglected (but the ec¬ 
centricity damping is included), the results match the observa¬ 
tions much better. However, no mechanism capable of suppress¬ 
ing type I migration over the whole inner disk has ever been 
found. Recent studies have shown that MRI-driven disk winds, 
in which gas material is blown away from the surface of the disk, 
can alter the density profile of the gas disk, potentially slow¬ 
ing down or even reversing the migration of the protoplanets 
(e.g.,[^zuki & Inutsukal2009HSuzuki et al.l2010tlC)gihara et all 

l2015h . This possibility, however, requires further investigation. 
Unless a mechanism for a global reduction of type I migration is 
demonstrated, our results imply that in situ accretion of close-in 
super Earth is unlikely. 
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